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Abstract 

Colombeau's generalized functions are used to adapt the distributional approach 
to singular hypersurfaces in general relativity with signature change. Equations 
governing the dynamics of singular hypersurface is obtained and it is shown that 
matching leads to de Sitter space for the Lorentzian region. The matching is possible 
for different sections of the de Sitter hyperboloid. A relation between the radius of 
S 4 , as the Euclidean manifold, and the cosmological constant leading to inflation 
after signature change is obtained. 
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1 Introduction 



Hartle and Hawkingfl], in constructing a satisfactory model of the universe, try to avoid 
the initial spacetime singularity predicted by the standard model of cosmology using a 
combination of general theory of relativity and quantum mechanics. The basic features 
of the so called " Hawking Universe" obtained are as follows: 

1. A satisfactory theory of quantum gravity will represent the gravitational field, in 
the manner of general theory of relativity, by a curved spacetime[l]. 

2. The proper understanding of ordinary quantum mechanics is provided by Feyn- 
man's "path-integral" or "sum-over-histories" interpretation. In ordinary quantum 
mechanics the basic idea is that a quantum particle does not follow a single "path" 
between two spacetime points, and so does not have a single "history", but rather 
we must consider all possible "paths" connecting these points. Therefore the usual 
wave function ^ is interpreted as an integral over all possible "paths" that a quan- 
tum system may take between two state. To solve the path integral, however, one 
must rotate the time variable in the usual quantum mechanical wave function to 
imaginary values in the complex plane, which yields as the new time coordinate the 
"Euclidean" time r = it [1,2]. 

3. There is a wave function for the entire universe that is given by a Feyn- 
man path integral. The basic idea here is that one sums over all possible four- 
dimensional spacetimes( or spacetime "histories") connecting two three-dimensional 
spaces (states). In order to evaluate the path integral, however, one must again ro- 
tate the time variable to imaginary values which changes the integral from Lorentzian 
to Euclidean one. The result is that the temporal variable in the wave function is 
changed to a spatial one. In other words, tyjj sums only over Euclidean spacetimes, 
that is, over four-dimensional spaces with positive definite signature (+ + + + 
)[l-3]. 

4. One wants to reach a certain state S at which the evolution of the universe becomes 
classical, in accordance with general theory of relativity and standard model of 
cosmology. Accordingly, Hawking proposes a path integral over the Euclidean four- 
space g^ u , and matter-field configurations that yields S. S is characterized by the 
three-metric /i^- and a value of the scalar field, <p. 

5. To avoid an initial spacetime singularity, the cosmic path integral will include only 
compact (or closed) four-geometries, so that the three-geometries, marking succes- 
sive states of the universe, shrink to zero in a smooth, regular way [6]. Hawking's 
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universal wave function is obtained, therefore, by integrating only over compact 
four-geometries (Euclidean "spacetimes" ) that have the 3-space S as the only (lower) 
boundary and are such that a universe in state S will subsequently evolve. 

Statements (3) and (5) above are the essence of the idea that the universe was initially 
Euclidean and then, by change of signature of "spacetime" metric, the transition to usual 
Lorentzian spacetime occured. Earlier attempts to describe this interesting aspect was 
based on Euclidean path integral formulation of quantum gravity and the analogy to 
quantum tunelling effect in quantum mechanics [1]- [6]. 

The rise of this idea led many authors to consider it within the classical theory of general 
relativity [7-20], with some controversies regarding the nature of the energy-momentum 
tensor of the hypersurface of signature change[15,18,19]. Ellis and his coworkers[8], have 
shown that classical Einstein field equations, suitably interpreted, allow a change of sig- 
nature of spacetime. They have also constructed specific examples of such changes in 
the case of Rabertson- Walker geometries. Ellis [9], by constructing a covariant formalism 
for signature changing manifolds, has shown the continuity of geodesies over the chang- 
ing surface. Hayward[10] gives the junction conditions necessary to match a region of 
Lorentzian-signature spacetime to a region of Riemannian-signature space across a space- 
like surface according to the vaccum Einstein or Einstein-Klein-Gordon equations. As we 
will show, Hay ward's junction condition R = is not compulsory. He has also considered 
the increasing entropy, large-scale isotropy, and approximate flatness of the universe in 
the context of signature change [11]. Dereli and Tucker [12], described classical models of 
gravitation interacting with scalar fields having signature changing solutions. Kossowski 
and Kriele[13] considered smooth and discontinuous signature change and derived neces- 
sary and sufficient junction conditions for both proposals. They investigate the extent to 
which these are equivalent. They have claimed that non-flat vaccum solutions of Einstein 
equations can only occur in the case of smooth signature change. Hellaby and Dray[15, 
19] argue for a non-vanishing energy-momentum tensor of the signature changing surface. 
Hayward disagrees to this result and favour a vanishing energy momentum tensor [18]. 
Kriele and Martin[16] do not accept the usual blief that signature change could be used 
to avoid space-time singularities, except one abandon the Einstein equations at the sig- 
nature changing surface. Martin[17] has studied Hamiltonian quantization of general 
relativity with the change of signature. He has also studied cosmological perturbations 
on a manifold addmiting signature change. 

We intend to use the distributional method of Mansouri and Khorrami [21,23] to approach 
the signature changing problem within the general theory of relativity. This is a powerfull 
formalism which can also be adapted to our problem, using the Colombeau's generalized 
theory of distributions [24-28]. This remedies the difficulties of non-linear operations of 
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distributions in the framework of classical theory of Schwartz_Sobolev[29]. In section 2 
we give an overview of the Colombeau's generalized theory of distributions. Sectoin 3 
contains the adaptation of the distributional formalism for singular surfaces in general 
relativity to the case of signature change, followed by concluding remarks in section 4. 
Conventions and definitions: We use the signature ( — h ++) for Lorentzian region and 
follow the curvature conventions of Misner, Thorne, and Wheeler (MTW). Square brack- 
ets, like [F], are used to indicate the jump of any quantity F at the signature changing 
hypersurface. As we are going to work with distributional valued tensors, there may be 
terms in a tensor quantity F proportional to some 5-function. These terms are denoted 
by F. 



2 A short review of Colombeau theory 

Classical theory of distributions, based on Schwartz-Sobolev theory of distributions, doesn't 
allow non-linear operations of distributions [29]. In Colombeau theory a mathematically 
consistent way of multiplying distributions is proposed. Colombeau's motivation is the 
inconsistency in multiplication and differentiation of distributions. Take, as it is given in 
the classical theory of distributions, 

e n = e v n = 2,3,..., (i) 

where 9 is the Heaviside step function. Differentation of (H) gives, 

n9 n - x 9' = 9'. (2) 



Taking n = 2 we obtain 
Multiplication by 9 gives, 
Using (2) it follows 



299' = 9'. (3) 
29 2 9' = 99'. (4) 

-9' = -9', (5) 

3 2' v ' 

which is unacceptable because of 9' ^ . The trouble arises at the origin being the 
unique singular point of 9 and 9' . If one accepts to consider 9 n ^ 9 for n = 2, 3, . . ., the 
inconsistency can be removed. The difference 9" — 9 , being infinitesimal, is the essence of 
Colombeau theory of generalized functions. Colombeau considers 9(t) as a function with 
"microscopic structure" at t — making 9 not to be a sharp step function ( Fig.l), but 
having a width t[24]. 9(t) can cross the normal axis at any value e where we have chosen 



4 



it to be e > §[29]. It is interesting to note that the behaviour of 9(t) n around t = is 
not the same as 6>(t)(Fig.2), i.e. 6(t) n ^ 9(t) around t = 0[24]. In the following we give a 
short formulation of Colombeau's theory. 

Suppose $ G D(R n ) with D(R n ) the space of smooth(i.e.C°°) C-valued test functions on 
R n with compact support and 

J <&(x)dx = 1. (6) 
For e > we define the rescaled function $ e (x) as 

= (7) 

Now, for / : R n — > C, not necessarily continuous, we define the smoothing process for 
/ as one of the convolutions 

/»:= J f(yMy-x)d n y, (8) 

or 

f t (x):= J f{yW{y-x)d n y. (9) 
According to (7), equation (9) has the following explicit form 

fe(x) := J f(y)U(^)d n y. (10) 

This smoothing procedure is valid for distributions too. Take the distribution R, then 
by smoothing of R we mean one of the two convolutions (8) or (9) with / replaced byi?. 
Remember that R is a C-valued functional such that 

$ G D(R n ) (R, $) G C, (11) 

where (R, $) is the convolution of R and $. 

Now we can perform the product Rf of the distribution R with the discontinuous func- 
tion / through the action of the product on a test function \1>. First we define the product 
of corresponding smoothed quantities R t with f t and then take the limit 

(Rf, *) = Jim J R e (x)f e (x)^(x)d n x. (12) 

The multiplication so defined does not coincide with the ordinary multiplication even 
for continuous functions. Colombeau's strategy to resolve this defliculty is as follows. 
Consider one-parameter families (f e ) of C°° functions used to construct the algebra 

E M {R n ) = {(/ e ) | f e G C°°(R n ) \/K C R n compact, 

Vet eN n 3N e N, 3r]>0, 3c> (13) 
such that sup \ D a f e (x)\ < ct~ N V0 < e < 77}, 
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where 

D a = (14) 

(dx 1 )^ ■ ■ ■ (dx n ) a - ' 1 ; 

and 

\a\ = a\ + a 2 + h a n . 

Accordingly, C°°-functions are embedded into ^(iR n ) as constant sequences. For con- 
tinuous functions and distributions we require a smoothing kernel (j)(x), such that 

J d n x (f(x) dx = 1 and J d n x x a (p(x) =0 |a| > 1. (15) 

Smoothing is defined as (10) for any function /. 

Now, we have to identify different embeddings of C°° functions. Take a suitable ideal 
Af(R n ) defined as 

N(R n ) = {(/ e ) | (/ e ) G S M {R n ) VK C M n compact, 

Va G N n , \/N G N 3n> 0, 3 c> 0, (16) 
such that sup xeK \ D a f t (x)\ < c <E N V0 < e < n}, 

containing negligible functions such as 

fix) - J d n y^( y -^)f(y). (17) 
Now, the Colombeau algebra Q(M n ) is defined as, 

««■) = ^ (18) 

A Colombeau generalized function is thus a moderate family (f e (x)) of C°° functions 
modulo negligible families. Two Colombeau objects (/ e ) and (g e ) are said to be associate 
(written as (g e ) ss (/ e )) if 

lim e ^o Jd n x (f e (x) -g £ (x)) p{x) =0 V y? G D(J? n ). (19) 

For example, if = </?(— x) then 56* ~ |<5, where 5 is Dirac delta function and is 
Heaviside Step function. Moreover, we have in this algebra 6 n 6 and not 9 n = 9. For 
an extensive introduction to Colombeau theory, see [24,25]. 
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3 Distributional Approach to Signature Change 



There are two methods of handling singular hypersufaces in general relativity. The mostly 
used method of Darmois-Israel, based on the Gauss-Kodazzi decomposition of space-time, 
is handicapped through the junction conditions which make the formalism unhandy. For 
our purposes the distributional approach of Mansouri and Khorrami (M-Kh) [23] is the 
most suitable one. In this formalism the whole space time, including the singular hyper- 
surface, is treated with a unified metric without bothering about the junction conditions 
along the hypersurface. These conditions are shown to be automatically fulfilled as part 
of the field equations. In the M-Kh-distributional approach one choose special coordi- 
nates which are continuous along the singular hypersuface to avoid non-linear operations 
of distributiuons. Here, using Colombeau algebra, which allows for non-linear operations 
of distributions, we generalize the M-Kh method to the special case of signature changing 
cosmological models. 

Consider a spacetime with the following FRW metric containing a steplike lapse function 
ds 2 = -f(t)dt 2 + a 2 (t)( 2 + r 2 d9 2 + r 2 sin 2 6dtp 2 ), (20) 

where 

f(t) = 6(t)-9(-t). (21) 

It describes a signature changing spacetime with the singular surface t — 0. The metric 
describes a Riemanian space for t < and a Lorenzian space-time for t > 0. As has been 
argued before fit), analogous to 0(t), has a microscopic structure around t = with a 
jump equal to r as shown in Fig. 3. We choose again 

9{t) |i =0 = e with e > -. (22) 

Since 9(-t) = 1 - 9(t), we have 9(-t) = 1 - e and 

/(t)=2e-l. (23) 

This value gives us the correct change of sign in going from t < to t > 0. This 
"regularization" of f(t) at t — allows us to use operations such as /(t) -1 , f{t) 2 , and 
|/(t)| -1 . The physical interpretation of this behaviour of f(t) is that the phenomenon of 
signature change occurs as a quantum mechanical tunneling effect. This tunneling occurs 
in a width equal to the width of the jump, i.e. r. 

In what follows we consider f(t) to be the regularized function f e , defined according to 
Colombeau's algebra. Now, we are prepared to calculate the dynamics of the signature 
changing hypersurface in the line of M-Kh procedure [23]. First we calculte the relevant 
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components of the Einstein tensor, G u and G rr , for the metric (20): 

3 -2a f + fa l f( 3 -2a f + fa 3 2a fa - a fa + 4(kf 2 + fa 2 ) } ( . 
Utt ~ 2. fn. + 2 J \2 f2n 2 n2f2 /' ^ 



2 fa T ^2 Pa 2 a 2 f 
and 

p _ 1 2aaf-aja+Afa 2 +Akf 2 1 a 2 f 3 -2a/+g/ 3 2aaf~afa+A{kf 2 +fa 1 ) \ /r>r\ 

Lr " 1 ~~ 2 / 2 (l-fcr 2 ) 2 1-fcr 2 \2 / 2 a 2 _R 2 / 2 / V 

According to standard calculus of distributions, we have 

/(*) = 0(*)-0(-t) 

= £(*) + £(-*) = 2 <$(*), 

taking into account 5(— f) = . Now, using Colombeau algebra we can write 



(26) 



B(t)6(t) « |*(*). (27) 



Therefore we may write 



/(f) S(t) = 9(t)5(t) - 9(-t)5(t) 

« W*) - |*(*) (28) 
« 0, 

In evaluating (24, 25) we should take care of the following property of association. Having 

AB « AC, 

we are not allowed to conclued 

B « C. 

Now, using the relations (27) and (28) we obtain for the singular parts of equations (24) 
and (25) 

G tt = 0, (29) 

^-j^j^m- (30) 

where multiplication of the distribution 5{t) with the discontinuous functions is defined 
as in (12). 

According to [23] the complete energy-momentum tensor can be written as 

Tijlv = 0(f) T+ + 9(-t) T~ + CS, U 5(t), (31) 

where T^ v are energy-momentum tensors corresponding to Euclidean and Lorentzian re- 
gions respectively and C is a constant which can be obtain by taking the following pill-box 
integration defining 5 MI/ [23]: 

= lim / (T^u - g^u-)dn = - lim / G^dn, (32) 

2^ — ^0 J — S K K ^ — *0 </ — S 
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Since 

% v = CS^mx)), (33) 

and 

J % v dn = CS, V J 5{§(x))dn = CS^\^\, (34) 

we find 

£7=1^1 = 1^^*1. (35) 

where <3> = t = defines the singular surface E. The vector n M is normal to the surface $ 
and n measure the distance along it. Using the metric (20) we obtain 

c = my (36) 

The distributional part of the Einstein equation reads as follows: 

Gfx V = kT^. (37) 
we obtain using equations (29,30,33,36,37): 

= j^jj Su6(t), (38) 

and 

Now using equation (12), we must define the multiplication of ^-distribution with the dis- 
continuous functions jjj and -p. To this end we consider them as Colombeau's regularized 
functions, 

G u {t) = 5 e (t)(jj^) (40) 

and 

G*(t) = 6 e (t)(j^) . (41) 

Now according to (12), these two multiplications are as follows, 

1 



(*(*) |^ > *) = J™ / Gu(tMt)dt (42) 



and 



m-^fy *) = Jim / G 2e(t) *(t)dt (43) 

for any test function, Now we argue that G± e and G 2e are associate in Colombeau's 
sense, i.e. 

lim / (G u (t) - G 2e (t))*(t)dt = 0. (44) 
9 



This is correct because although G± e and G*2e are divergent at a common point, the dif- 
ference in their " microscopic structure" at that point tends to zero for e — > 0. Therefore, 
we obtain from (38), (39), and (44) the final form of the energy- momentum tensor of the 
singular surface, or the dynamics of, S: 

S tt = 0, (45) 

and 

S ~ = -4T^)- < 46 > 

If, along the above line, we compute the See and S </,</, we find, 

S„ = (47) 

K 

and 

K 

Taking the metric (20) and the results (45-48) we obtain for the so-called 'energy-momentum' 
tensor of the singular hypersurface 

kS£ = diag(0, -2H, -2H, -2H) , (49) 

where, as usual, we have chosen 

H = -. (50) 

CI 

Integrating this equation, gives 

a = a e Ht , (51) 

which shows de Sitter expansion of scale factor corresponding to the inflationary phase 
after signature change. This is in accordance with the features of the so-called "Hawking 
Universe" [1]. 



4 Junction conditions and interpretation of the re- 
sults 

Take the metric (20) with a defined as in (50, 51): 

ds 2 = Tdt 2 + e 2Ht (dx 2 + dy 2 + dz 2 ). (52) 

These coordinates cover half of the de Sitter Hyperboloid for — oo < t < +oo(30). The 
hypersurfaces t — to are spacelike, except for t = —oo which is null and the boundary of 
coordinate patch. 
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The Euclidean section can be interpreted as a S* 4 with radius R = In the 

Lorentzian section we are faced with an exponentially expanding universe in which the 
Hubble parameter is the constant H — ^, irrespective of signature changing time t . The 
cosmological constant in the Lorentzian sector is A = 3H 2 , where Therefore, there is the 
following relation between the cosmological constant and the radius of S 4 : 

A = | (53) 

Now, the boundary of the Euclidean part is defined by T = constant. These sections are 
S* 3 having a radius equal to a = exp HT. If the matching is along a section corresponding 
to less (more) than a hemisphere then the radius a of S 3 will increase (decrease) with T [1]. 
Hence, the expansion in Lorenzian part, i.e. a > 0, means that the matching corresponds 
to a part of S 4 which is less than a hemisphere where the radius of the boundary increases 
with the coordinate T. 

We now go on to calculate the extrinsic curvatures to check the junction conditions [23]. 
In the Lorentzian region the non-vanishing components of affine connection are 

r?i = r° 2 = r° 3 = He 2H \ (54) 

The extrinsic curvature is defined as 

K tj = efe^rv, (55) 

where e*, the mutualy normal unit 4- vectors in signature changing surface $, are defined 

as 

e t - g i , i - 

are coordinates adopted to the signature changing surface E. V M shows the covariant 
derivative with respect to the 4-geometry. We then find for the non-vanishing component 
of extrinsic curvature 

K~ = K 2 " 2 = K~ = He 2H \ (56) 

or 

K' 1 = K 2 2 = Kz 3 = H fort>-oo, (57) 

and 

K' 1 = K 2 2 = K^ 3 = fort = -oo. (58) 
Similarly, if we consider plus sign in (52), we find for the Euclidean embedding 

K+ 1 = K+ 2 = K+ 3 = -H, (59) 
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independent of t being < t < it. The jump in the extrinsic curvature is thus 

[K{] = K' j - Kt j = 2H, fort>-oo, (60) 

and 

[K{\ = H, for t = -oo. (61) 
Now, the junction condition for the matching is [23] 

KSl = \Ki]-H\K] (62) 

Comparing (60,61) we see that the matching condition is satisfied only for t = — oo. We 
therefore conclude that for these coordinates chosen for the de Sitter space the matching 
to a Euclidean space is possible only for the section of space representing t = — oo, which 
is a null hypersurface. 



5 Matching in a different coordinates 

The previous section shows that the requirement of signature change leads us to the de 
Sitter space in the Lorenzian part of the manifold and S 4 in the Euclidean part. This is in 
agreement with the results of quantum cosmology and minisuperspace considerations [1]. 
It is however well known that different coordinates for de Sitter space corresponds to 
different t = constant sections. The above de Sitter coordinates, for example, which 
cover only half of the de Sitter hyperboloid, are equivalent to k = flat three space 
sections[30] . 

Now, there is another useful coordinates familiar from Robertson- Walker metrics which 
cover the entire de Sitter hyperboloid and its t = constant sections are surfaces of constant 
curvature k — 1. These S 3 spaces have a minimum radius equal to We will now try 
the matching along these sections of de Sitter space. Consider the following metric with 
approprtiate lapse function corresponding to signature change [30] 

ds 2 = -f(t)dt 2 + a 2 cosh 2 (a-H) (d X 2 + sin 2 X (d9 2 + sin 2 6d4> 2 )) (63) 

where f(t) is defined as in (21). Following the same procedure as for (20) and again 
using Colombeau's algebra, we fined for elements of energy-momentum tensor of the 
hypersurface 

kS v = diag(0, tanh(a _1 t), tanh(a _1 t), tanh(a _1 t)), (64) 

a a a 

The hypersurface of signature change is defined as t = t . The no n- vanishing compopnents 

of the extrinsic curvature are 

K~ i = -itanh((rt ) i = 1,2,3. (65) 
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The corresponding components in the Euclidean embedding are 




tanh(a 1 T ), 



(66) 



where we have used a different symbol T as the coordinate in the Euclidean region corre- 
sponding to 'time'. This is allowed because the hypersurface of signature change as the 
boundary of the two different manifolds is defined differently in each of them. Now the 
jump of extrinsic curvature on the signature change surface is 



Now the junction condition (62) is satisfied for t = T = 0. Therefore, the matching is 
done along the section of de Sitter space with the minimum radius y^, which is equal to 
the radius of corresponding section of the Euclidean manifold. 

6 conclusion 

As we have seen in the previous sections signature change is only possible for de Sitter 
space in the Lorenzian part of the manifold. Depending of the coordinate patches one use 
on the de Sitter hyperboloid, different matchings are possible. Using de Sitter coordinates 
which cover only half of the hyperboloid the matching is along the t = — oo which is a 
null hypersyrface. For the Robert son- Walker type metric covering the entire hyperboloid 
the matching is along the hypersurface of minimum radius of the de Sitter hyperboloid. 
In both cases the cosmological constant is given by the radius of the Euclidean S 4 . We 
have therefore a geometrical interpretation for the value of the cosmological constant. 
Vanishing of the cosmological constant would mean that the Euclidean manifold is an R 4 
which is noncompact and therefore contradicts the quantum cosmological considerartions. 
Now, it is very reasonable to take for the radius of S* 4 the only fundamental length that 
we have, namely the Planck length. We therefore come to the conclusion that Ipi = Jjr- 
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